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^ Abstract. We investigate the structure of the nonequiUbrium stationary state 

1/^ , (NESS) of a system of first and second class particles, as well as vacancies 

(holes), on L sites of a one-dimensional lattice in contact with first class parti- 
cle reservoirs at the boundary sites; these particles can enter at site 1, when it 
is vacant, with rate a, and exit from site L with rate /3. Second class particles 
can neither enter nor leave the system, so the boundaries are semi-permeable. 

el The internal dynamics are described by the usual totally asymmetric exclusion 

' process (TASEP) with second class particles. An exact solution of the NESS 

was found by Arita. Here we describe two consequences of the fact that the 
j^ , flux of second class particles is zero. First, there exist (pinned and unpinned) 

fat shocks which determine the general structure of the phase diagram and of 
the local measures; the latter describe the microscopic structure of the system 
at different macroscopic points (in the limit L — > oo) in terms of superpo- 
sitions of extremal measures of the infinite system. Second, the distribution 
of second class particles is given by an equilibrium ensemble in fixed volume, 
or equivalently but more simply by a pressure ensemble, in which the pair 
potential between neighboring particles grows logarithmically with distance. 
S ^ We also point out an unexpected feature in the microscopic structure of the 

w . NESS for finite L: if there are n second class particles in the system then the 

distribution of first class particles (respectively holes) on the first (respectively 
last) n sites is exchangeable. 
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1. Introduction 
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. . In recent work Arita [U [2], using a matrix ansatz, found the nonequilibrium 

C^^ ' stationary state (NESS) of a new version of the widely studied one-dimensional 

totally asymmetric exclusion process (TASEP) [3]-[TT] (see in particular 12J for a 
recent review of matrix methods for the TASEP) . The model is defined on a subset 
of the one dimensional lattice Z consisting of L sites. Each site i, i = 1,...,L, 
may be occupied by a first class particle, occupied by a second class particle, or 
vacant; vacant sites are also referred to as holes, and first class particles simply as 
^j ' particles. We shall let these three possible states correspond to the values 1,2, and 

0, respectively, of a random variable r^; we also introduce the indicator random 
variables r]a{i), a = 0, 1, 2, such that r]a{i) = 1 if r^ = a and ria{i) = otherwise. 

The internal (bulk) dynamics of the system are given by the usual rules for 
the TASEP with second class particles ^. The occupation variable r^ at site 
i,i = 1,...,L— 1, attempts when r^ = 1 or 2 to exchange at rate 1 with Ti+i; 
when Ti = I the exchange succeeds iff Ti+i = or 2, while for r^ = 2 it only 
succeeds if tj+i = 0. In other words, a first class particle at site i jumps to the 
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right by exchanging with either a hole or second class particle at site i + 1, while 
a second class particle can only jump if the site on its right is empty. At site 
i = 1, first class particles enter the system at rate a provided that site 1 is vacant 
(ri = 0); at site i = L, first class particles leave the system at rate /3 provided 
that site L is occupied by a first class particle (r^ — 1). Second class particles 
are thus trapped inside the system; since only first class particles can cross the 
boundaries, we refer to these as semi-permeable. (A similar system, but with a 
different form of semipermeable boundary, was considered in [I3|). An equivalent 
system is obtained by interchanging first class particles with holes, left with right, 
and a with f3, and this symmetry will be reflected in the structure of the NESS. 
The latter will be determined by the parameters a, (3 and the density 7 = n/L of 
second class particles, where n is the number of second class particles in the system. 
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Figure 1. The cross section of the phase diagram at a fixed 7. 



The phase diagram of this system in the limit L — > cx) is given in Figure [T] [2] . 
The diagram is determined by the distinct formulas for the (first class) particle 
current Ji in the different regions: 



r 1-7' 



(1.1) 



for a, (3 > ac (region I), 



Ji = < 



a(l — a), for a < ac, ol < [3 (region II), 



, /3(1-/?), for /? < ac,/3 < a (region III); 
here the critical value a^. of a and (3 is 

(1-2) "^^^- 

The current J2 of the (trapped) second class particles must vanish and the current 
Jo of holes satisfies Jo = — Ji. Note that the form of Ji as a function of the 
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parameters characterizes the phase plane regions, and that Ji is continuous but 
not smooth across all the boundaries. 

The phase diagram is similar to that of the open one component TASEP [71 [H] 
and indeed in the limit 7 — > reduces to it; moreover, in regions II and III the 
current is independent of 7 and takes the same values as in the one component case, 
although the size of these regions shrinks as 7 increases. As we discuss in Section|6l 
however, there will be residual differences between the local microscopic states of 
the one species model and 7 — > limit of the two species model; in particular, there 
remain an infinite number of second class particles near one or both boundaries of 
the two species system. Note also that there is a discontinuity, equal to 7, in the 
derivative of Ji with respect to a (/3) on the I/II (I/III) boundary; one might say 
that the order of the phase transition in Ji when 7 7^ differs from that when 
7 = 0. 

The macroscopic density profiles Pa{x) in the NESS, a = 0, 1, 2, defined by 

(1.3) Pa{x)^ lim (?7aW), 0<a;<l, 

L — >oo, n/ L — >7, i/ L — >x 

with (•) the expected value in the NESS, have been computed in 2 ; the results 

are summarized in Table [1] (but see Remark II . 1 1 below) . Knowing any two of these 

densities determines the third, via X]a/^a(^) ~ ^- ^'^ fact, knowing pi{x) for all 

X and all values of a and j3 determines po{x), from the particle- hole symmetry, 

and hence all profiles, but for clarity we give in Table [T] both pi{x) and pq{x). In 

region II the system divides itself into two parts, x < xq and x > xq, with different 

formulas for p(){x)] similarly in region III there are different formulas for pi{x) for 

X < xi and x > xi. Here 

7 
xo = 1 - -; r-, a < /3, a < ttc; 

(1.4) ^'-'" 

xi = i_2B ' P-^^ P<ac. 

On the II/III boundary a = (3 < ac, the shock line, the profiles include linear 
regions: 

(1 — a) H a, < X < xq. 



(1.5) 



Pq{x) ^ \ xq xo 

Xii <x <\ 

a, < X < xi, 

Pl{x) = { l-X X-Xl //I 

■aH (1 — a), xi < x < 1. 



1 — xi 1 — xi 

These arise from averaging over the position of a shock, as in the one species TASEP; 
further discussion is given below. 

Remark 1.1. (a) The density values Pa{x) at the boundaries a; = 0, 1 and at the 
fixed shocks a; = xq, xi may depend the way the limit (|1.3p is taken. The boundary 
cases x = 0, 1 were discussed in [2] except on the I/II and I/III boundaries. We 
discuss the limits at xo,a;i in Section [5] and [51 this gives some further information 
about limits at the boundaries since xq = and xi = 1 on the I/II and I/III 
boundaries, respectively. 

(b) In the one-component model the phase plane regions corresponding to I, II, 
and III are called the maximum current, low density, and high density regions. 
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Table 1. Density profiles in different regions of tlie phase plane. 
Note that xq is defined only in region II and on its boundaries, and 
xi only in region III and on its boundaries. 



Region 


Pi{x) 


Po{x) 


I 

I/II boundary 

I/III boundary 


OLc 


ac 
ac 

ac 




X <X\ 


X > Xl 


X < Xo 


X > Xl^ 


II 

III 

II/III boundary (Shock Line) 


a 


1-a 


a 


P 


1-/3 


P 


a{=P) 


linear 


linear 


a {=(3) 



respectively. We do not adopt that terminology here since in region III the particle 
density is low for x < xi. 

Some typical profiles, obtained from simulations, are shown in Figure [2] Since 
these are from a finite system, they do not coincide perfectly with the description 
in Table [TJ there are boundary effects, and the density transitions in regions II and 
III, at Xl and xq respectively, have nonzero width of order yL. This is related to 
the nonuniqucness of the limit (jl.3p at these points, as mentioned above, and is 
discussed in Section [5l 

We now give an intuitive discussion of some of the phenomena that give rise to 
these profiles. Consider a macroscopically uniform portion of our system in the 
limit L — > cxD, with densities of holes, particles and second class particles denoted 
by po I Pi and p2 , respectively, where po + Pi + P2 = 1 • In such a region the measure 
will be the known translation invariant measure, with these densities, for the two 
species TASEP (see [9l [Ml [T^ and the discussion in Section [8]) ; in this measure 
(which is not a product measure) the first class particles considered separately, and 
the holes considered separately, are distributed according to product measures, so 
that Ji — pi{l ~ pi) and Jq — — po(l — Po)- Thus from Jq = —Ji it follows that in 
any uniform stretch of the NESS either 



(1.6) 



Pi =Po = (1 -P2)/2 



P2 = 0, Pi = 1 - Po- 



This fact, which may be seen in the results of |2J, is key to understanding the gross 
structure of the densities in different regions of the phase diagram. 

These density profiles differ from those the single-species open TASEP in two 
notable ways: in regions II and III the density profiles have a point of discontinuity, 
and on the II/III boundary the linear region occupies only part of the system. These 
and other properties can be understood in terms of the occurrence of a fat shock. 
By this term we refer, not to a broadening of the sharp shock usually seen in the 
TASEP (as can occur in the partially asymmetric model [16]), but rather to a 
macroscopically uniform interval which contains all the second-class particles, and 
thus conforms to the first alternative in (|1.6p . We may think of this fat shock as 
bounded by shocks of the usual sort occurring in two different one species TASEP 
systems which are obtained by making appropriate identifications of second class 
particles with either first class particles or holes. To see how this occurs, recall that if 
one either identifies first and second class particles by coloring holes black and both 
kinds of particles white, or else identifies second-class particles and holes by coloring 
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REGION II: a = 0.30; (1 = 0.50; T= 0.25 



REGION III; (1 = 0.50; p = 0.20; y= 0.26 
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I / II /III TRIPLE POINT; a = 0.375; p = 0.375; 7= 0.25 
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II /III SHOCK LINE; a = 0.30; P = 0.30; y= 0.25 



Holes; Second class; First class; 



0.25 x„ 0,50 X, 0.75 

(d) 

REGION I; ct= 0.50; p = 0.50; 7= 0.25 



0.25 0,50 0.75 1 

(e) 
Figure 2. Density profiles in a system with L — 1000. 




these species red and first-class particles blue, then the black/white particles, as 
well as the red/blue particles, form standard two species TASEPs in the bulk. The 
dynamics at the boundaries is different, since some white "particles" or red "holes" 
will be trapped in the system. A careful justification of the conclusions below is 
given in Sections [3] and [5] 

Consider now the behavior of the system on the boundary of regions II and III 
(the shock line). Then by previous analysis, see e.g. [TD], one knows that a typical 
profile for the one species model contains a shock between a region of density a on 
the left and 1 — a on the right; the shock position has mean velocity zero and its 
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(fluctuating) position is uniformly distributed over the system. We see this same 
behavior for both the black/white and red/blue systems described above, with the 
black/white shock necessarily located to the left of the red/blue one. The typical 
profile at any given time looks on the macroscopic scale like Figure [3l where the 
convention is that at any point x the height of the region labeled with particle type 
a is paix). 



Holes (black, red) 



Second-class 
particles 
(white, red) 



1-a 



Particles (white, blue) 



Figure 3. Shock interpretation at a = f3 < ac- Densities Pa{x) 
are plotted against x, with the convention that at the height at x 
of the region labeled with particle type a is pa {x) . The fat shock 
may in fact be located anywhere in the system. 



Clearly both shock fronts have mean velocity zero and are trapped in the system, 
and since the total number of second-class particles is 7L, the macroscopic width 
w of the fat shock must satisfy w = 7/(1 — 2a). This forces the two shock fronts to 
move (i.e., fluctuate) in collusion so as to keep the macroscopic shock width fixed; 
we expect this fluctuation, as for the shock in the single component TASEP on its 
shock line, to be on a diffusive time scale growing as L^. The density profiles Pa{x) 
arise as averages over the shock position, and this gives rise to the linear profiles 
(jl.Sp : in contrast to the situation in the one species case, however, here they occupy 
only part of the system because the shock can fluctuate only over an interval of 
width 1 — w. The shock fluctuation is also reflected in the structure of the local 
measures obtained in the limit L ^ 00, which are superpositions of states with 
different densities (see Section[5|). The critical value of a occurs when the fat shock 
fills the system, i.e., when w = 1, from which we regain (jl.2p . 

The situation in regions II and III is similar. The fat shock width is in general 



(1.7) 



w{a,P,j) = 



1 



l-2aA/3' 

where a /\ (3 — min{a,/3}; in region II the shock is pinned to the right boundary, 
and in region III is pinned to the left boundary. Since the shock is fixed it gives 
rise to discontinuities in the density profiles; see Figure [2] as well as the discussion 
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of a related model, where similar behavior occurs, in Section [8] (Figured]). There 
is no corresponding discontinuity in the single-species TASEP (7 — 0) because in 
that case there is a single shock of zero macroscopic width. In region I, w as given 
in (|1.7p is greater than 1 and the fat shock fills the system; the density profiles are 
uniform and conform to the first alternative of (jl.6p . 

The outline of the rest of the paper is as follows. In Section [2] we discuss the 
matrix method for this system. We use a different representation of the matrices 
from that of [5], which makes it easier to prove certain features of the NESS dis- 
cussed later. In Section [3] we show that the marginal distribution induced by the 
NESS on particles in the first n states of the system, and on holes in the last n, is 
exchangeable, i.e., that the probability of finding r first class particles (holes) on 
some specified set of r sites among the first (last) n is independent of the choice of 
sites. In Section |4] we establish the fat shock picture described above. In Section [5] 
determine the local measures, in the bulk, for the infinite volume limit of the sys- 
tem, and in Section [6] consider the local measures near the boundaries, focusing on 
a Bernoulli property which is a consequence of the exchangeability established in 
Section [3l In Section [7] we show that the second class particles form an equilibrium 
system, most simply described by a pressure ensemble. This is related, in our case, 
to the fact that the current of second class particles is zero. For similar situations, 
see [nillH]. 

In Section [8] we make some concluding remarks and, in particular, describe some 
closely related models. One such model is a generalization of the standard "defect 
particle" model; another describes a system of first class particles, second class 
particles, and holes on a ring with one semi-permeable bond which second class 
particles cannot cross. Several more technical remarks are recorded in the appen- 
dices. 



2. The matrix ansatz 

The stochastic system described in Section [T] is ergodic in finite volume L and 
thus there exists a unique invariant measure /i^'^ on the configuration space 

(2.1) YL,n = { (n, . . . , tl) I Tj = 0, 1, 2; n ^ 2 for n values of i }, 

where from now on we will assume that < n < L. This measure may be obtained 
from a matrix ansatz [2\ , combining the matrix algebra of [9J (which discussed the 
system with the same constituents as in the current work, but on a ring) with 
the treatment of the one species open system via matrix-elements from j7]. One 
introduces matrices Xa, Xi, and X2 and vectors \Vp) and {Wa\ which satisfy 

(2.2) XiXq = Xi + Xq, ^1^2 — X2, X2XQ — X2, 
and 

(2.3) {W^\Xo = -{W^\, X^\Vf3)^hvp). 

a p 

Then for a configuration t — [n^- ■ ■ ,tl) (z YL,n the probability of r in the invariant 
measure is 

(2.4) {T)..,=Z'^-''mn)-\W^\Xr,---Xr,\Vf3), 
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where Z°''^{L,n) is the normaUzation factor 

(2.5) Z"'0{L,n)^ Y. {Wo,\Xr,---Xr,\Vp), 

which, with a shght misuse of the nomenclature of equihbrium statistical mechanics, 
we call the partition function. We will frequently omit superscripts such as a, (3 in 
(j2.5p when no confusion can arise. 

We will work throughout in a realization of (|2.2|) - (|2.3p . different from that of 
[2j, for which the matrices and vectors have the further properties 

(2.6) X2 = XiXo-XoXi = \Vi){Wi\, X2\Vp)^\Vi), {W^\X2 = {Wi\, 

(2.7) (Wa^lVi) = {Wi\Vp) = 1 for all a,/3, 

where \Vi) and (W^il are vectors satisfying (|2.3p . Note that X2 is then a one- 
dimensional projection operator. The realization is given in Appendix \^ but we 
will need no consequences beyond (|2.6p and (|2.7p . Because of (|2.7p we make the 
convention that ^"'^(0,0) = Z'^''^{0,0) = 1. 

Remark 2.1. The nature of X2 in this representation shows that certain distribu- 
tions obtained from (|2.4p factorize. Let Qi, . . . , Qn denote the (random) positions 
of the second class particles in the system (note that these can be ordered once and 
for all). Then the probability that the jf^ , j^, . . . , j^ second class particles are 
located on sites qj^, . . . , qj^^ is 

f^tniQj. =qj^,t^l,...,m)^ Z^^P{L,n)-'Z^-\q,, - 1, ji - 1) 

m 

(2.8) X J] Z^^\q,^ - q,^_, - 1, J, - j,_i - 1) Z^^^{L - q,^ , n - j„). 

i=2 

Moreover, if we condition on the event that Qj = g, i.e., that the j'*^ second 
class particle is located at site q, then the conditional measure is a product of the 
measures associated with the sites before and after j, so that if t is a configuration 
consistent with this event then 

(2-9) f^t'.ni'^l, ■■■,rL\Qj =q) = fiq\j^liTl, . . . ,T,_l)/i^'^^„_^.(T,+ l, . . .,Tl). 

A factorization property of this type is also known for the translation invariant 
measures for the two species TASEP [9j [15] . Similar expressions are easily obtained 
when conditioning on the presence of several specified second class particles at 
specified sites. We will use ()2.8p and ()2.9p in Sections [4] and [5l when we discuss the 
fat shock and describe local measures in the NESS. 

3. Exchangeability of Measures 

In this section we demonstrate a remarkable property of the finite- volume NESS 
with n second class particles: the exchangeability [12] of the measure on first class 
particles within the first n sites, or equivalently on holes in the last n sites. Specif- 
ically, this means that for any r < n the probability of finding first class particles 
on the r sites 1 < ii < i2 < ■ ■ ■ < ir < n depends only on r, i.e., is independent 
of the choice of positions ii,i2, ■ ■■ ,ir- When r = 1 this is implicit in (38) of ^2], 
although it is not emphasized. As a consequence of the ideas of the proof we will 
also obtain, for any i,j with i,j > 1 and i + j — 1 < L, the probability of finding 
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a block of j consecutive first class particles starting at site i; this generalizes the 
density formula of [2], which corresponds to j = 1. 

The key quantity for our arguments is the probability of finding (first class) par- 
ticles at sites ii, . . . ,ir-i together with a block of j particles starting at site ir, 
where r > 1, ii < • • • < i,_i, and v > ir-i + 1; we allow j = 0, with the inter- 
pretation that in this case there is no restriction on what happens at site v or any 
succeeding sites. Thus the probability in question is Z(L, n)~^Er{L, n; ii, . . . , ir]j), 
where Er is a sum of weights for certain configurations t G Y^n'- 

(3.1) Er{L,n-ii,...,ir;j)^ ^ {Wo.\Xr, . . . Xr,\Vp) 

ri^=---=ri^_^=\ if r > 2 
T-ir=n,, + l=---=-ri,,+j-l = l if j > 1 

In p.ip we must have v+j~l < i, since there are only L sites, and r+j — 1 < L — n^ 
since there are n second class particles. For certain parts of the analysis we will 
have to consider separately two cases, in which these two inequalities respectively 
provide the effective bounds on j: 
Case 1: v > " + r, so that < j < L — i,. + 1; 
Case 2: i,. < n + r — 1, so that {) < j < L — n — r + 1. 
We will analyze the Er using a simple recursion: 

(3.2) Er{L,n;ii, . . . , v;i) = Er{L,n;ii, . . . ,ir: j + l)+Er{L-l,n;ii, . . . , v; j-1). 

This holds whenever all terms are defined, which requires that j be positive and 
satisfy j < L — ir in Case 1 and j < L — n — r in Case 2. To derive p.2p we consider 
the value of ti^+j in each term of the sum (|3.ip . Terms with t^+j — 1 sum precisely 
to EriL,n;ii,.. .,ir-; j + 1), ar 
matrix algebra to reduce 

the resulting sum is just Er{L — l,n;ii, . . . ,ir;j — 1). 

To determine the Er the recursion ()3.2p must be supplemented by boundary 
conditions at the maximum and minimum values of j. When j = 0, p.ip gives 

fo A\ T? fT ■ ■ n\ \Er-i{L,n]ii,...,ir-i;l), ifr>2, 

(3.4) Er{L,n;ii,...,ir;Q) = < 

Z(L, n), II r = 1. 



terms with r^-i-j 


= or Ti+j = 


2 we use the 


1^0 = Xi + Xq, 


if Ti^+j == 0, 




1^2 = X2, 


ifn,+, =2; 





The value of Er for j maximal is case dependent. In Case 1, if j takes its maximum 
possible value L — ir + l then each matrix product in p.ip ends with X-^ '''^^iVg) = 

^-(L-v-fi)|y^^^ SO that 

Er{L,n;ii, . . . ,ir;L - ir + 1) 

. . _ rr(^-''■+^)i?,.-l(v-l,n;^l,...,v_l;l), if r > 2, 

^^■^^ -\p-(^-^^+^^Z(^r~l,n), ifr = l, ^^^'' ^^^ 

In Case 2, if j takes its maximum possible value L — n — r + 1 then the rightmost 
factor in the matrix product in (|3.ip is X2 and there are no factors of Xq; using the 
matrix algebra relations X1X2 — X2, X^ = X2 and X2|Vg) — |Vi), we have that 

(3.6) Er{L,n;ii, . . . ,ir; L — n — r + 1) = 1 (Case 2). 
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Lemma 3.1. The recursion ()3.2|) . together with the boundary conditions (13. 4|) and 
either p.5p or (|3.6p . determines Er{L,n;ii, . . . ,ir',j) by an inductive computation. 

Proof. The primary induction is on increasing values of r, with n held fixed through- 
out. The inductive assumption that the Er-i are known is needed when (|3.4p or 
p.Sp is applied, and since for r = 1 the right hand side of these equations is a 
(known) partition function one may treat all values r > 1 uniformly. Then for 
fixed r (r > 1) and ii, . . . , v we induce on increasing values of L: for the minimum 
possible value, L = i^, we must be in Case 1 and have j = or j = 1, so that 
all Er{L, n;ii, . . . , ir]j) are determined by (|3.4p and (|3.5p . For any larger value of 
L either (j3.5p or (|3.6p determines Er{L,n;ii, . . . ,ir',j) for the maximum possible 
value of j and we may then induce downward on j using ()3.2p . D 



We can now verify exchangeability; we show that Er{L, n;ii, . . . , ir]j) is equal 
to the corresponding weight with the sites ii, . . . , v in standard positions 1, . . . , r. 

Theorem 3.2. Let 1 < ii < 12 < ■ ■ ■ < ir be sites such that ir < n + r — 1, and let 

j be an integer less than or equal to L — n — r + 1. Then 

(3.7) Er{L,n;ii,...,ir;j) = Ei{L,n;l;r + j - I). 



Proof. The proof is by induction on r. By Lemma 13.11 it suffices to show that 
-Ei(L,n;l;r + i — 1) satisfies the same recurrence relation (J3.2I) and boundary 
conditions (|3.4p . (|3.6p as does Er{L, n;ii, . . . , ir\j)- This follows immediately from 
the corresponding relations for Ei and, for r > 2, the induction assumption. D 

We finally give the explicit formula for Ei which, by (|3.7p . also provides a formula 
for Er{L, n;ii, . . . , ir] j) when v < n + r — 1. This result will not be needed in the 
remainder of the paper. 

The formula involves the Catalan triangle numbers [20] 

(3.8) ^^^fm + n\n.~n+l 



n J m + 1 
which satisfy the recursion 

(3.9) c™-i+c:L, = cr 

and the boundary conditions 

(3.10) c!!\ = o, cj' = i, c;;: = c™_i = - 

We then define the additional constants 

(3.11) c,.=(f.' '[^:^'- 



2m 
m 



m. 



(3.12) dmjk = <^ " ' ' ""^ 

1^4,0, if j = m + 1. 

In (|3.1ip - (|3.12p the convention is that (^) = for integer p, q with p > 0, g < 0. 
Theorem 3.3. Case 1: If n + 1 < i < L and < j < L + l-i then 

L — i L — i—1 

(3.13) Ei{L,n;i;j) = Y,Cj.kZ{L~k,n) + Z{i~l,n) ^ dL-»j,fc/5"^+'+''-^ 

k=j fc=0 
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Case 2: If I < i < n and < j < L — n then 

L — n 

(3.14) Ei{L,n;i;j)= J^ ''o,kZ{L - k,n), 

k=j 

Proof. Case 1: We temporarily denote the right hand side of (|3.13p by F{L, n; i; j). 
By Lemma 13.11 it suffices to verify that F satisfies relations corresponding to 
()3.4p . (|3.5p . and (|3.2p . It will be convenient to denote the two terms in p.l3p 
by Fi{L,n;i;j) and F2{L,n]i\j), respectively. 

Since dL-i,o,k = and co,fc = S^fl we have immediately that F{L, n; i; 0) = 
Z{L, n) (compare p.4p ). Moreover, the sum defining Fi{L, n; i; L — i + 1) is empty 
and so from d^-i^L-i+i^k — ^k.o we have F{L, n; i; L — i + I) ~ /3^^^'^^^Z{i — 1, n) 
(compare p.Sp ). It remains to check the equivalent of (13.21) . which we shall show 
is satisfied by Fi and F2 separately; recall that in (|3.2p . 1 < j < L — i. For j > 2, 

Fi(L,n;z;j + l) + Fi(L-l,?i;i;j-l) 

L—i L—l — i 

= E CtU^{L-k,n)+ J2 CtUZ{L-l-k,n) 

k='j + l k=j-l 

L-i 

(3.15) =E (^fe-j-i +^M) Z{L-k,n)^F,{L,n;i;j), 

k=j 

where we have used Cti^ = (see ([3Jn| 'l and C^l}_i + C^Z] = Clz] (see dSU). 
The case j = 1 is easily checked separately. Similarly one verifies that 

(3.16) F2{L,n;i;j + l)+F2{L-l,n;i;j-l)^F2{L,n;i;j) 
separately for j < L ~ i -- 1 and for j = L — i. 

Case 2: We denote the right hand side of (|3.14p by G{L,n;i;j), and show that 
G satisfies the appropriate boundary conditions and recursion. One checks im- 
mediately that G{L,n;i;0) = Z{L,n) (compare p.4p ) and, using (j3.10p . that 
G{L, n;i;L — n) = 1 (compare p.6p ). Finally one shows that, for 1 < j <i — n— 1, 

(3.17) G{L,n;t;j + l) + G{L-l,n;t;j~l)^G{L,n;i;j); 

the proof is essentially the same as that of the recursion for Fi in Case 1. D 

4. The fat shock 

In this section we give a precise definition and analysis of the fat shock discussed 
informally in the introduction. The analysis will be used in the next section for the 
determination of local states in the infinite volume limit. We define the fat shock 
microscopically as the region between the positions Qi and Q„ of the first and last 
second class particles in the system. 

The joint distribution of Qi and Qn was obtained in Remark 1 2. 11 it is convenient 
here to write this, for j, fc, Z > 0, as 

02;f (j, k, I) = filfXQl = J + l.Qn=J+k + 2) 5,+k+LL-2 

(4-1) Z'^'^mn) '5,+.+,i_2. 

We can determine the large-L behavior of 9 by replacing the partition functions in 
(|4.ip with their asymptotic forms; these can be obtained from [7] and [2], and are 
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sumniarized in Appendix [BJ In some cases it is convenient to further approximate 
the distribution of k, which represents the fat shock width on the microscopic scale, 
by a Gaussian. (RecaU that a macroscopic width w = w{a,f3,^) for the fat shock 
was predicted on heuristic grounds in Section 1 (see (|1.7p ). so we expect that k '^ wL 
for large L.) We omit details of the computations and summarize the results in the 
next remark; where the notation 9'^'^{j,k,l) ^ f{a,/3,j,j,k,l) indicates that the 
ratio of the two quantities approaches 1 as L — > cx) with n — [7LJ , where [uj is 
the greatest integer contained in u. 

Remark 4.1. (a) On the boundary of regions II and III {a — P < ac), 



(4.2) e^Jj,k,l)^j^^^—^^-^e > ' S,+,+i,L-2, 

where A{a) = (1 — 2a)'^/(47a(l— a)). That is, under O^'nU^ ^' 0: ^ is approximately 
Gaussian with mean Lw and variance of order L, j is approximately uniformly 
distributed on the range < j < L — I — 2, and I — L — 2 — j ~ k. On the 
macroscopic scale, this means that the width of the fat shock is w and its left 
endpoint is uniformly distributed on the interval [0, 1 — lu]. 

(b) In region II {a < a^a < l3), 

(4.3) etiij, fc, I) ^ p-^^~-)'^{l) yMe-^(")('=-^-)V^ <5,+.+,,_2, 
where we have introduced the (normalized) probability distribution 

(4.4) pi^il)^ — uZ^ (i,0), / = 0, 1,..., 

defined for u < /3(1 — /3) if /? < 1/2, u < 1/4 otherwise. p'"'^{l) decreases expo- 
nentially for large I unless u — 1/4, when the decrease is as l~^/^ (see (jB.!!) ): p is 
normalized by (|B.2[) . Thus on the microscopic scale I is typically of order 1 and the 
shock width k is distributed as in (a). On the macroscopic scale the fat shock has 
width w and is pinned to the right end of the system. The analysis in region III is 
similar. 

(c) On the boundary of regions I and II (o^c = a < /3), 



+k+l,L-2- 



(4.5) eliU, fc, /) ^ p"(i-")-'3(0 2\/^e-^(")('=-^)'/^ 5, 

' V ttL 

This is as in (b) except that here w — 1 and as a result k is distributed as a Gaussian 
random variable conditioned to have value at most equal to its mean, and there is a 
corresponding factor of 2 in the normalization. The analysis on the I/III boundary 
is similar. 

(d) At the triple point {ac = a ~ P), 

(4.6) etiU, k, I) ^ j^^^ e'^i'^K'^-^f/L s^^,^,,_,. 

The distribution of k is as in (c) but here j and I are free, subject only to the 
constraint j-\-l = L — 2— k. 

(e) In region I {ac < a, (3), 

(4.7) 0li{j, k, I) ^ p(i-^^)/4."(j) p(i-^^)/4./3(0 5,+,+,^_2; 
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j and / are both of order 1 (microscopically) and k is determined by the delta 
function constraint. 

Note that the results of Remark 14.11 confirm the picture of the fat shock behavior 
sketched in Section 1. 

5. Local states in the infinite volume limit in the bulk 

In this section we discuss a question inspired by the treatment of the one com- 
ponent system by Liggett (lOj : is there a local state at position x of the system 
(in the infinite volume limit), and if so what is it? To formulate a precise question 
we consider a limit in which n and i increase with L in such a way that i — > cxd, 
L — i ^ cxD, i/ L — > a; G [0, 1], and n/ L — > 7 g (0, 1). In this setting we ask about 
the existence and nature of the weak limit Voul^oo 2^~*Ml nj where T is translation 
by one lattice site and so the operator T"' carries site i to the origin; equivalently, 
we consider the sites of our open system to run from \ — i to L — i and look at the 
probabilities of configurations in the interval from —K to K, take L, i, L — i, and n 
to infinity as above, and then make K arbitrary. The limit (if it exists) is a measure 
on the configuration space Y — {0, 1,2}^; we call it a local state in the bulk since 
(in the L ~* 00 limit) it describes a situation infinitely far from each boundary; the 
local state at the boundary will be discussed in Section [6] 

It will suffice to consider a special class of these limiting procedures; specifically, 
we will always take 

(5.1) n = iiL = Lt^J ^-^.d i = iL = \xL\ + cvi; 

we must assume that c > if x = and c < if a; = 1. We then define 

(5.2) M.,c^ lim r-^>2;^^. 

The limit in (j5.2p certainly exists along subsequences, by the compactness of the 
set of measures on Y . To simplify notation we will ignore the necessity of passing 
to subsequences; since the limiting measure will be found to be unique, the limit 
of the sequence itself must also exist. For most values of the parameters the limit 
will in fact be independent of the choice of c, but this is not true when a; = xo in 
region II or on the I/II boundary, or a; = xi in region III or on the I/III boundary. 
We first consider the currents and densities in the state ^x.c- The currents in 
the finite system, and hence also in the limit, are independent of the site: 

(5.3) (ryi(j)(l - 771 (j + 1)))^^^ = ((1 - r^^{j - l))%(j))^,., = Ji 

for any j, with Ji given in (jl.ip . The limiting densities Pa{x, c) for a = 0, 1, 2 are 
defined by 

(5.4) pa(a;,c) = lim (77a(u)) „,0 =(77^(0)) 

with the last equality expressing the fact that i^ corresponds to the origin in fj,x,c- 
It is easy to check, from the asymptotic computations of [2] , that the limit in (|5.4p 
would be unchanged if iz, were replaced by II + j for any fixed j, which implies 
that (rjaij)) n = Pa{x,c) for any j, i.e., the densities under fix,c are translation 
invariant. Equation (|5.4|) may be viewed as a refined version of (|1.3p . in which the 
ambiguity in the L -^ 00 limit there has been removed. 

Noting that in the L —f 00 limit the generator of the dynamics in the neighbor- 
hood of iL does not involve any boundary terms or any constraints on the densities 
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of the three species beyond J^aVaU) = 1, one verifies easily [10] that Hx.c must 
be invariant for the infinite-volume two species TASEP dynamics. It then follows 
that fix,c must be a convex combination of the extremal invariant measures for the 
infinite volume two species TASEP. These measures have been classified in 14]: 
there is (i) a family of translation invariant measures j/'^f^-^i^ defined for Aq, Ai > 0, 
Ao + Ai < 1, in which holes, first class particles, and second class particles have 
densities Aq, Ai, and 1 — Aq — Ai, respectively, and (ii) a family of non-translation- 
invariant "blocking" measures j)™'"^ where ttt,, n G Z U {— oo, cx^}, m < n, and m, n 
are not both infinite: j>™'" is a unit point mass on the configuration r™'" given by 

0, if i < m, 
(5.5) t"'" = <! 2, iim<i<n, 

ii n < i. 

However, the translation invariance of the densities implies that none of the blocking 
measures can be present in the superposition giving ^^^c- 

Thus there exists a probability measure K^.ddXo, dXi) (which depends also on a, 
/3, and 7) that specifies the weights of the different translation invariant extremal 
measures which enter into the superposition: 



(5.6) fJ^x,c= / Ka;,c(dAo,dAi)i/ 

JAo,Ai>0, Ao-I-Ai<1 

We will see that: (i) for most values of the parameters, Kx,c is a point mass, so that 
fj,x,c is one of the measures j/-^0''^i , (ii) in some cases, in which x may lie to the left 
of, within, or to the right of the fat shock discussed in the introduction, fix.c is a 
superposition of the two or three measures corresponding to these possibilities, and 
(iii) no more complicated superposition can occur. Note that, as a consequence, 
the current Ji is the same for all elements of the superposition (and the same holds 
for Jo and for J2 = 0). Here is a first result in this direction. 

Theorem 5.1. If fJ,x,c is defined by (|5.2p and the current and densities at x are re- 
lated by Ji = po{x,C){l-po{x,c)) = pi{x,c){l-pi{x,c)), thenpx.c = i/Al(2;,c),pi(a:x)_ 

The condition in the theorem that po(a;, c)(l — Pa{x, c)) = Pi{x, c)(l — pi(x, c)) 
corresponds to the zero current of second class particles and leads to the alternatives 
of p.6p . We see from Table[T]that this theorem determines Px,c completely for most 
but not all values of a, j3, 7, x, and c and that the results are consistent with the 
intuitive picture sketched in the introduction. In summary: 



Remark 5.2. It follows from Theorem 15. II that: 

(a) In region I the local state px,c is i/"-'^^""-, in particular, it is independent of x and 



c. 

,1 — a,Q 



(b) In region II the local state px,c is zy^~"'" for x < xq and v°'-°' for x > xq, i.e., 
respectively outside and inside the fat shock. Region III is similar: the local state 
is u^'^ for x < xi and ijl^^^~l^ for x > xi. 

Other cases are not determined by the theorem: 

(c) The local state is not determined by Theorem 15.11 in the interior of the regions 
where either type 2 or type particles have a linear profile, that is, on the II/III 
boundary (where the fat shock is not pinned to one or the other end of the system) 
with < X < xo or 1 > a; > xi. See Figure [2l[c,d). 
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(d) The local state is not determined by Theoreni l5.ll (i) at a; = xq in region II and 
on the I/II boundary, where xq — 0; (ii) at x = xi in region III and on the I/III 
boundary, where xi = 1, or (iii) at a; = xq = and a; = xi = 1 at the triple point. 
All of these points are edges of the (pinned) fat shock; see Figure [2l[a,b,e). 

We will determine ^x,c in cases (c) and (d) below. 

ProofofTheorem l5.lt Using ()5.6p together with the relations (??i(«))j,ao,ai — Ai, 
{r]i{i){l — ?7i(i + l))),yAo,Ai = Ai(l — Ai) (which hold for all i because the marginal 
of iy^o,>^i on configurations of first class particles is a product measure), we find 
that 



(5.7) pi(x, c) = (??i(«))^^_^ ^ r"^ KxAd^o,d\i) = (Ai)^^_^ , 

and 

Ji = (r;i(z)(l-,7i(* + l)))^^_^ 

(5.8) = /"Ai(l-Ai)K,^,(dAo,dAi)-(Ai(l-Ai)),^^. 

From Ji — pi{x,c){l — pi{x,c)), then, we see that (A^) = (Ai)^ , so that 
Ai = (Ai)^ = pi{x,c) almost surely with respect to k^.c- Similarly, Aq = po{x,c) 
almost surely with respect to n^.c, so that p^.c — i/P°'^^''^^'P^''^''^\ D 

We now turn to the determination of the local measure px,c at a point x where the 
densities are varying linearly or are discontinuous — case (c) or (d) of Remark 15.21 
Recall that in Section [4] we have determined the probability d^n^j, k, I) that Qi = 
j + 1 and Qn = j + k + 2, where j + k + l = L — 2 and Qi and <5„ are the position 
of the first and last second class particles. Now let Ml'„ ,■ k i denote the measure 
f^Ln conditioned on Qi = j + 1, Qn = j + k + 2. The key observation we will use 
follows from Remark 1 2. 1[ specifically, from (|2.9p or a simple generalization thereof: 
if G is a function on Yl^u which depends on the Ti only for toq < i < mi, then 



•(G),=^, iimi<j, 

(5.9) {G)^^.0 . . = i ^^'''^'^) t^l-\ , ' if J' + 2 < "lo, ?7ii < 2 - 

{T-J-^-^G) I,,, ifj + fc + 3<mo. 



Now let {tl) be a sequence of integers with r^ ^ oo and tl/ \L — > 0. For any 
function F on Y depending only on finitely many spins we may write 

(F)^ = lim y eli {j,k,l)(T'^F) ^,, 

j+k+l=L-2 ' 

(5.10) = ^hni^ [S« (F)^a) + sf (F)^,., + S^^^^ (F)^<3, + remainder] . 

Here //]^ is for p = 1, 2, 3 the probability measure defined by 



(5.11) p^l^^E^^^-'J2''letiAj,k,i)T'^^pZ 
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where J^) 1 1 ranges over values satisfying j > il + fL, ^2) L i o^*^^ J < ih — r^ and 

,(3) 



k > iL + tl, J2) k I °^^^ k <iL - tl, and 



(5.12) 



^(P) 



^(p) = y^^^^ Q^-P (n h n - ) ,,"./5 



V"l('3i > *i' + ''i)' ifp=l, 



j,k,l 



;l;,(j,fc,0 = <M:;j:(Qi<«L-^L, Q„>iL + rL), ifp = 2, 



tJ-nliQ" < ^L-rL), ifp = 3. 



The remainder in (|5.1ip contains those terms from the fuU sum over i and k which 
are omitted from ^ , ^ , and X) • We have suppressed the dependence of 
these various entities on a, /3, 7, x, and c. 

We now take the L ^ 00 hmit in (IS.lOp . It follows from Remark 14.11 and 
the fact that r^ grows more slowly than yZ that the remainder vanishes in this 
limit. The S]^ are expressed as probabilities in (|5.12p and their limiting values 

Si?c = liniL^oo 2)^ may be determined from Remark 14.11 these limits will be 
discussed on a case by case basis below. 

Finally, to study lim^^oo Ml '^6 observe that for sufBciently large L (if F de- 
pends on Ti only for |i| < m then r^ > m suffices) we have by (15. 9p that 

'4'^"'ES,.^Ll(j,fc,0(T-F)^.., ifl>=l, 

(5.13) (F)^^(,) = hf~'i:flieliSj,k^){T^--^-^F)^,^^^_^^, ifp=2, 

[^f'" ES, oil 0' fc' (r--^-^-^F),.^ , if P = 3. 

The limits limL-»oo /^^ ^^^ -P ~ l7 2,3 are all treated similarly; let us discuss the 
case p = 2 in detail. Equation ()5.13p displays ^)^ as a convex combination of the 
probability measures 7'~^*^~"'~"'^Vfe'n-2- Each of these measures, for large L, is by 
Remark 15. 2r a) approximately equal to u'^^l^^°'^P (recall that a A /? = min{a,/3}), 
since the critical value a* of a for a system with k ~ wL = 7^/(1 — 2a A /3) sites and 
n ~ 7L second class particles — and thus an effective value 7* = n/k — 1 — 2a f\ (3 
of 7 — is (1 — 7*)/2 = a A /3. The same should be true oi ii\ . The corresponding 
evaluations for p = 1,3 come from the results of |10] for the local measures in the 
one species open system. We conclude that 

{j^l-QA(l/2),aA(l/2)^ ifp=l, 
j^aA/3.aA/3^ if p = 2, 

^/3A(l/2),l-/3A(l/2)^ ifp = 3, 

and so 

(5.15) ^ix c = S(.l)iyl-"^(V2),aA(l/2) ^ 2(^2)^^aA/3,QA/3 ^ ^^(3) ^/3A(1/2)4-;3a(1/2) _ 

Remark 5.3. The heuristic argument for (|5.14p given above could be made precise 
by justifying the implicit exchange of limits; we sketch instead an alternate proof, 
again for p = 2. We know that the limiting current for the measures /ifc'„_2j as 
L and hence k ~ wL goes to infinity, is a(l — a), and the limiting densities of 
holes, particles, and second class particles are a, a, and 1 — 2a, respectively. One 
can in fact show further that these limits are obtained with error which goes to 
zero uniformly at sites i satisfying r^ < i < k — tl, from this, it follows that the 
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(2) 

measures /i}^ have the same hmiting current and densities. Then an argument as 
in the proof of Theorem 15.11 estabhshes (|5.14p . 

To complete our discussion of the local states in the bulk we must supplement 
()5.15p with a determination of the weights Ex,c = limi^oo ^l ^^^ cases (c) and (d) 
of Remark 15.21 The cases in the next remark are parallel to those of Remark 14.11 

Remark 5.4. (a) On the boundary of regions II and III (the shock line, case (c) 
of Remark l5.2p we find from Remark l4.ir a) and (|5.13p that 

({l-w~x)+, ifp=l, 

(5.16) EiPl = —^ X h-w~{l-w-x)+-{x-w)+, ifp = 2, 

[{x-w)+, ifp = 3. 

Here u+ — u ii u > and u_|_ = if w < 0. Note that these coefficients, and hence 
the local measure fix,c given by (|5.15[) . are independent of c. 

(b) In region II [a < ac,a < /3), at the fixed shock xq, the SS do depend on c: 



5(3) =0. 



Xo,C 



(5.17) SW, = 1 - $ (cv/aR) , 42), = $ (cv/IR) , 
Here $ is the error function defined by 

(5.18) $(i) = ^ / e-"'/2^r. 

The analysis in region HI is similar. 

(c) On the boundary of regions I and II (ac = a < /3), where xq = 0, we find that 
for c > 0, 

(5.19) 4^^ = 2-2$(cv/IR), 4'}^2^[cy/A(a))-l, sg^ = 0. 
The analysis on the I/HI boundary is similar: for c < 0, 

(5.20) e['} = 0, e['1 ^ 1-2^ [c^/A(a)), S^^^ = 2$ (c^AR) . 

(d) At the triple point (qc = ct ~ P), where a;o = and xi = 1, (|5.19p and (|5.20p 
again hold (with c > and c < respectively). 

We finally note that the results of this section yield density profiles as well as the 
finite volume corrections to these profiles at the fixed shocks (see Remark 1 1.1 f a)), 
since pa{x,c) may be calculated from (|5.4p and (|5.15p . For example, on the II/III 
border (shock line) we find in this way that 

(5.21) po(^)-sW(l-«) + (sE + 4'^)«, Pi(a:) = (sW+S(> + 43)(l~a), 



with the weights S^p^ given by (|5.16p : it is easy to see that (|5.2ip reproduces (|1.5p . 
Here we have used the notation Pa{x) of (|1.3p . rather than writing pa{x,c) as in 
(|5.4p , since the densities do not depend on c. In region II we have, from (|5.15p and 
dnUl), that 

(5.22) po(a;o,c)-l-a-$(cv/I(^) (l-2a). 

Similar results hold in region HI at the point xi. 
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6. Local states in the infinite volume limit at the boundaries 

In this section we study limiting measures limL^oo T '^Ml « ^^ in (|5.2p . still 
taking n^ — L7^J but now assuming that cither i^ or L ~ i^ is fixed. Without loss 
of generality we can assume that i^ = 1 or i^ = L (the measure as seen from site 
J or site L — j can be recovered from these limits) and thus define 

(6.1) Wcft= lim T^Vif^Li, Mright = lim T^-^/x^f .,. 

L — >oo ■ L ' -I L — ^oo ■ L ' -I 

Note that /iicft (respectively /iright) does not coincide with any of the measures /io.c, 
c > 0, (respectively /xi^c, c < 0,) studied in Section[5l The densities under /iioft and 
Mright were studied in [2]; for example, (?7a(i))u ■ is denoted pj^^jj ■ in 0. 

By the particle hole symmetry it suffices to consider /iicft , which is a measure on 
the semi-infinite configuration space {0, 1, 2}^^'^'^'^''^^--}. In general we do not have a 
proof that the limit defining fiicu exists (except along subsequences), although we 
expect this to be true; see also the comment below Theorem 16.11 The next result, 
however, gives a somewhat surprising property which any (subsequence) limit must 
satisfy; to simplify notation, we will speak as if the limit itself exists. 

Theorem 6.1. The distribution of first class particles under the measure /iieft is 
Bernoulli with a constant density p. where p is given by 



(6.2) 



Proof. By Theorem 13.21 we know that the (marginal) distribution of the variables 
771 (i) under /iioft is exchangeable, so that by de Finetti's theorem [19j this distri- 
bution is a superposition of Bernoulli distributions. From [^ we know that for 
any i > 0, p = (?7i («))„, ,. = hmL^oo ('7i(*))„°-3 is given by (|0)) and that 
limi^oo (77i(i)(l — rii(i + 1)) c,/3 = Ji (see (11.11) '). In each region of the phase 

plane these limits satisfy the relation J = p(l — p). Then from an argument as 
in the proof of Theorem 15.11 it follows that the ?7i(i) distribution is the product 
measure i''', where p is given by (|6.2p . D 

Note that in region II the density of second class particles any finite distance 
from the left boundary goes to zero as L ^ 00 [2], so that knowing that the 
distribution of particles is Bernoulli completely determines any limiting measure 
to be a Bernoulli measure on particles and holes only. It follows that the limiting 
measure exists without passing to subsequences. 

We discuss briefly one aspect of the measure /iicft in the limit 7 — > (note that 
we are taking this limit after the L — > 00 limit). Consider first region I; from 
Remark |4. If e) we see that the position Qi of the first second class particle in the 
system is distributed according to p^/^'"((?i); this is a normalizable distribution 
which decreases as q^ for large qi, so that there remains a second class particle 
in the system, but ((9i)„ = 00. In fact more is true; by a calculation similar 
to that of Remark 14.11 one can show that all Qj — Qj-i, j = 2, 3, , . . ., have this 
same distribution (see also the discussion of the pressure ensemble in Section [7]) so 
that there remain an infinite number of second class particles in the system under 
Zoloft. The same is true in Region III, but there by Remark l4.ir b) Qi is distributed 



"c, 


in region I, 


a. 


in region II, 


/3, 


in region III. 
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according to p/3(i-/3),a^ g^ ^j^g^^ ('3i)u,„ft < °o; the distribution of the Qj — Qj-i, 
j — 2,3, , . . ., is the same as in Region I. 

Remark 6.2. One may compare Theorem 16.11 with resuh in 9 for the infinite 
volume Umit of a two-component TASEP system on a ring: that the distribution of 
first class particles to the right of a second class particle, and that of holes to the 
left of such a particle, is Bernoulli. The two results are closely related, because if 
we set a = /3 = 1 in the open system then the matrix element {Wi \Xti ■ ■ ■ Xr^ \Vi) 
giving the weight of the configuration ti , . . . , t^ is [5] exactly the weight of the 
configuration 2,ti, . . . ,tl on a ring. Because the numbers of first class particles 
and of holes on the ring is fixed, and these numbers fluctuate in the open system, 
this does not establish an exact equivalence of the a = j3 = 1 case of Theorem 16.11 
to the result of [|9j; nevertheless, it is clear that the former is in some sense a 
generalization of the latter to values of a and (3 other than 1. But the result of [3] 
is in another sense more general than Theorem 16. 1[ since the infinite volume limit 
of the open system has zero current of second class particles, but this is not true 
for the system of [9j. 

7. The pressure ensemble for second class particles 

We here consider the steady state distribution of the second class particles only, 
so that one may think of identifying the first class particles and holes as a new type 
of hole. For d a positive integer we define 

(7.1) 0„(d) = -log(4''^Z"'i(d- 1,0)) = -\og{A-''{WMX^ + X,Y-'\V,)) 

It follows from the {a, (3) symmetry of Z°'^^{L,n) that (t)p{d) is also equal to 
- \og{A-'^Z^'l^{d-l, 0)). Using (HU we find that the probability that the n second 
class particles in the systems are located at sites qi < q2 < ■ ■ ■ < qn ^^ 



where we have denoted <j>i{d) by (pid). The motivation for the factors 4 '^ in ()7.1|) 
will be discussed below; with this normalization (l){d) ~ (3/2) log d for d —i- oo 

mm- 

We note that (|7.2p has precisely the form of the canonical distribution for a 
system in a domain of length L with particles interacting with their nearest neighbor 
only via a pair potential 4>{d). (Such an interaction is rather unphysical; one may 
think of any intervening particle as screening the interaction of particles that it 
separates.). There is also a potential 4>a{d) {4>[3{d)) representing the interaction of 
the first (last) particle with the left (right) boundary. 

The TASEP dynamics for the full system gives rise in a natural way to a dynamics 
on the system of the second class particles which satisfies detailed balance with 
respect to this equilibrium measure. In the state in which the second class particles 
are at gi , . . . , (7„ the i"^ second class particle moves to site g^ + 1 at rate 1 whenever 
that site is empty (in the original sense), an event which by a simple generalization 
of (|2.9p occurs in the NESS with probability 



, . (M/i|Xo(Xo + Xi)^-+i-^--^|^i) ^ e-0(-?.+i-g.-i) . 
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and with probability 

One finds similarly that the probability that the site g^ — 1 is occupied by a first 
class particle is 

g-0(9i-?i-l-l) £-■01,(91-1) 

(7.5) —. r— , ifz>l, , I \ 1 ifi = l. 

The dynamics in which g^ ^ g^ + 1 when qi+i — Qi > 2, with rate given by (j7.3p . 
and gi — > gi — 1 when qi — Qi-i > 2, with rate given by (j7.5p . is easily seen to satisfy 
detailed balance with respect to the measure (|7.2p . 

To obtain the properties of the system described by (|7.2[) in the thermodynamic 
limit, L ~^ oo, n/L —^ 7, it is most convenient to consider the pressure or isobaric 
ensemble tt"'^ PUE^, in which instead of keeping the volume L of the system fixed 
we imagine that the right wall is in contact with a reservoir of pressure p. The value 
of p can be chosen so as to make the average volume equal to L, as discussed below. 
More precisely, we let the position of the right boundary, which we denote Qn+i, 
fluctuate, and add a term involving the pressure p to the measure. This yields the 
probability distribution in the pressure ensemble: 

TTp.'f (9l,---,9n,gn+l) 

(7.6) / " \ 
= Z'^-f\p, ny^ exp I -(l)a{qi) ~ ^ 0(gj - qi-i) - <pp{L - g„) - pq„+i j . 

The partition function has the form Z°''f^{p,n) = Zi{a,p)Z2{p)'^Zi{P,p), where 
Zi and Z2 are readily found, for z — y/l — e^P satisfying 

(7.7) 1 > z > max{0, 1 - 2a, 1 - 2/3}, 
to be given by 

7.8 Z^{a,p) = ^- '-, Z^ip) - -—. 

z + 2a— 1 1 + z 

Thus (|7.6p becomes 

(7.9) 
7r^;f (gi, . . . ,g„,g„+i) = Z,{a,p)~'e-^-^''^^-P'>^ 



Yl^2ipr'e~ 



{qi~qi-i)~p{qi-qi-i) 

.i=2 



2i{P,Pr e 



lp-</'^j(9" + l-?")-p('3..+ l-«n) 



The convenient factorization property of the probability tt" !^(qi, . . . , g„, qn+i) dis- 
played in (|7.9p . which implies that the variables qi and qj — qj-i, j = 2, . . . , n + 1, 
are independent, has prompted the use of the pressure ensemble for equilibrium 
systems, without any reference to dynamics. The requirement that particles only 
interact with their first neighbors is usually imposed artificially (see, however, [13]). 
In our model this condition arises naturally from the dynamics. Note that qn — Qi, 
the width of the fat shock, is thus represented as a sum of independent random 
variables; this is consistent with its Gaussian distribution in regions I and II of the 
fixed volume ensemble (see Remark H?T|) . 
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One easily checks that (writing now tt^ '^ = tt) 



dp ^ '"^ ' z{z + 2a- 1)' 
(7.10) ((jj-grj_i)^ = -— logZ2 = -, j = 2,...,n, 



dp " ''^^^ z(z + 2/3-1)' 

Note that when z approaches its lower bound, which is if a, /3 > 1/2 and max{l — 
2a, 1 — 2/3} otherwise, the average size (qn+i)^ of the system goes to infinity for 
every n > 1; there is simply not enough pressure to confine the system. To agree 
with standard definitions we have defined the potentials (f>a in ()7.ip so the the size 
of the system in the absence of boundary terms, that is, ((?„ — <?i)^, goes to infinity 
when p — > (z — > 0). 

To find the appropriate pressure corresponding to the canonical ensemble with 
L = n/j studied above we must set the expected system length 

(7.11) (g„+i), = (qi), + J2 (lj - <ij-i). + (9n+i - Qn)^ = -—logZ'^-^ip, n), 

i=2 ^ 

equal to L and solve for p (or z), subject to the restrictions ()7.7p . With (|7.10p the 
equation to be solved becomes 

a{l + z) n I3{l + z) ^ n 

^' ' z{z + 2a-l) z z(z-t-2/3-l) 7' 

We will discuss the solution of this equation in various regions of the phase plane; it 
is useful to bear in mind that each of the three terms on the left hand side increases 
as z decreases from 1 to its lower limit max{0, 1 — 2a, 1 — 2/3}. 

Consider first region I, where a, /3 > (1 — j)/2. Since for z = 7 the left hand side 
of (|7.12p is n/7 + 0(l), where the 0(1) term is positive, the solution must be of the 
form 2 = 7 + 0(l/n). In the limit n -^ 00 we thus have z = 7 or p = — log(l — 7^). 
In this case (gi)^ and all {qj — g^-i)^, j — 2, . . . ,n + 1, are of order unity, so the 
bulk of the system has, in the limit n ^ 00, the same structure as that obtained 
from our NESS when L — > cxd in region I. 

In Region II, where a < (3 and a < (1 — 7)/2, we have from z>l — 2a>l — 2/3 
that the third term in (|7.12p is 0(1), and from z > 1 — 2a > 7 that the first term 
must be 0(n), i.e., we must have z = 1 — 2a + 0{\/n). In fact we find easily that 
for large n, 

(7.3, ,.i_2a+Mi^i^ 

1 — 2a — 7 n 

Now (gi)^ is of order n but (g^+i — gi)^ is still of order 1; this corresponds to the 
fat shock being fixed to the right wall, i.e., to what we see in region II. The situation 
in region III is of course obtained by interchanging a with /3 and left with right. 
In the case when a = /3 < (1 — 27)/2 one sees again that z ~ 1 — 2a + c/n and 
that {qij^ = {q-n+i — (Inj^^ are both of order n; this simply means that the average 
position of the fat shock is in the middle. 

One may of course analyze the pressure ensemble directly, rather than looking 
for the correspondence with the open system of fixed length; the key question is 
how one allows p or equivalently z to vary with the number n of (second class) 
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particles. If z is held fixed (necessarily in the range (I7.7|l ) then the behavior of the 
system corresponds to that of the open system in region I. If a < 1/2 and a < /3, 
and one takes z = 1 — 2a + cjn then the behavior is as in region II; similarly if 
/3 < 1/2 and (] < a one obtains region III behavior by taking z = 1 — 2/3 + c/n, and 
if a = /3 < 1/2 such a z value gives behavior corresponding to the II/III boundary. 
A detailed analysis of the ensemble, for example of the shape of the profiles of the 
fat shock, would essentially repeat the analysis in the fixed L, i.e., fixed overall 
density 7, ensemble studied earlier, and we will not take up these questions again 
here. 

Note that for almost all permissible values of the pressure our system is in region 
I; only by fine tuning the pressure to change with n in a range of order 1/n do 
we get configurations as in regions II or III. This is reminiscent of what happens 
when one goes from a fixed magnetization to a fixed external magnetic field h 
in the Ising model at low temperatures, in dimension two or higher. The whole 
coexistence region, corresponding to the average magnetization being smaller than 
the spontaneous magnetization, corresponds to the single value h = 0. 

Other choices of z can lead to regimes different from those considered in the 
present work. For example, if we again suppose that a ~ (3 < 1/2, but now take 
z closer than order 1/n to 1 — 2a — to be specific, say z = 1 — 2q: + c/n^ — then 
((j„+i)^ is of order n^ and hence the density of second class particles is zero. 

8. Concluding Remarks 

1. As noted already in several places above, the local properties of our system away 
from the boundaries approach, as L ^ 00, those of the states of the two species 
TASEP on the lattice Z. Because of this it will be useful to describe here some 
known properties of the (extremal, translation invariant) NESS's of that system, 
i.e., of the states i/P^-p^ introduced in Section [5l These states differ from those of 
other models for which the NESS of the finite open system can be solved exactly, 
such as the one species simple exclusion process and the zero range processes [9l 
[MlEl], in that they are not product measures; this is so despite the fact that their 
projections (marginals) on the configurations of first class particles alone, or on 
the configurations of holes alone, are in fact Bernoulli. The states i/P°'P^ may be 
obtained [^ as the A^ — > 00 limits of states of a two component TASEP on a ring 
of N sites, with Na = PaN particles of type a, a = 0, 1, 2, where P2 = 1 ^ Po ^ Pi] 
see also [l4] . 

As noted in [9], the structure of the states vP°^p^ is quite intricate, containing 
several mysterious features which we still do not understand in any direct intuitive 
way. They are not even Gibbs measures [Tl], even though all the truncated cor- 
relation functions involving a finite number of sites decay exponentially fast. This 
decay follows from the (mysterious) fact that, conditioned on the presence of a 
second class particle at site i, the measure factorizes: the left and right sides of i 
become independent. For the corresponding property for the open system studied 
in this paper see Remark 12.11 Another (mysterious) fact is that if one conditions 
on the presence of a second class particle at i then the particles to the right of i, 
and the holes to the left of i, have this Bernoulli property [9]. 

Another related property of the states vP°'P^ is that if one conditions on there 
being a first class particle (respectively a hole) at site i then the measure to the left 
(respectively right) of i is the same as if there was no conditioning at all, i.e., the 
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same as that described in the first paragraph of tliis section. (This may be expressed 
cohoquially by saying that if one observes that the fastest horse is in front then one 
gains no information about the rest.) The property has in fact been established in 
not only the two species but also the n-species TASEP (see p6l, Proposition 6.2), 
using a representation of the stationary measure based on queuing theory; a direct 
proof for the two species model may be given using the two properties of second 
class particles noted in the previous paragraph. We remark that the property of 
factorization around a second class particle does not extend in a direct way to the 
n-species model ^^. 

2. The fact that the measures i/PO'P^ are not product measures gives extra struc- 
ture to the local states fix,c discussed in Section[5l which are superpositions of such 
measures. We note here however that as in the case of the one component TASEP, 
when such a superposition occurs only on the shock line a = (3 < 1/2, the transla- 
tion invariant measures which enter into the superposition (and which correspond 
to the measures on one side or another of a shock) all have the same current. This 
can be interpreted as saying that the properly averaged local current does not fluc- 
tuate. These averages can be obtained either as long time averages of the stochastic 
flux across a single bond, or as spatial averages over an interval of length K, with 
K ^ oo after L ^ oo. We believe that the convergence of the average total flux 
across an open system to a deterministic value, as L — > oo, is a general property 
of the NESS of systems like those discussed here, but do not know how to prove 
this directly at the present time. It seems reasonable to expect similar behavior 
in higher dimensions and different settings, e.g., for driven diffusive systems on a 
torus [27] . 

3. It follows from the "separating" property of conditioning on the presence of 
a second class particle at a given site that the distribution under i/P^'P^ of the 
second class particles alone is given by a renewal process [HI [H] ■ When the current 
J2 vanishes, i.e., when po = pi = [1 — P2)/2, then (as noted in Section [7]) the 
distribution of the distance between nearest neighbor particles in this process has 
a simple exponential dependence on p2 which can be obtained from a pressure 
ensemble, with p = — log(l — P2), as in region I of the open system. 

Combining this expression for the pressure as a function of the density with 
standard thermodynamic relations we can obtain expressions for the chemical po- 
tential A and Helmholtz free energy a in the uniform infinite system of second class 
particles with density p2'- 

A(p2) = -log' 



(8.1) VI + P2, 

a{p2) = (1 - P2) log(l - P2) + (1 + P2) log(l + P2)- 



From (j8.1[) we can obtain the large deviation function for the probability of finding 
n2C particles in an interval oi £ lattice sites [28] . The large deviation for first class 
particles or holes alone is of course given by the properties of the Bernoulli measure. 
Large deviation properties of the full measure have not, so far as we know, been 
determined for the two species system. 

4. Even knowing fully the properties of the infinite system still leaves open the 
problem of how fast the correlations in the vicinity of a site [_xL\ approach those in 
the local measure iix,c- This may be of particular interest in the case when ij,x,c is a 
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superposition of extremal infinite volume measures i/P°^p^ as discussed in Section [S] 
We might expect the L dependence in that case, where typical density profiles differ 
from average ones, to be different from that where the two coincide. We leave this 
as an open problem. 

5. We now briefly describe two related model systems, containing both first and 
second class particles on a ring, which are intermediate between those studied in 
[9j and in this paper. 

5.1 Recall the "truck" or "defect particle" model [Ml [301 [SI [12] , a standard two 
species TASEP system: a single defect particle together with (first class) particles 
and holes, on a ring of L + 1 sites, can exchange with a hole ahead of it (clockwise) 
at rate a and with a particle behind it at rate /3. Let us add to the ring also n 
(standard) second class particles, which make the same exchanges as does the defect 
particle but at rate 1 in each case, and which do not exchange at all with the defect 
particle. To be definite let us say that there are ni first class particles and no holes 
on the ring, with n + ni + no — L. Then the stationary measure for this system 
is almost the same as that of our open system: using the matrices Xq, Xi, and 
X2 of Section[2l and introducing also Xs = \Vf3){Wa\, we find that a configuration 
(5, Ti, . . . , Ti, where 6 represents the defect particle, has weight: 

(8.2) TviXsXr, ■■■XrJ= {Wo^lXr, ■ ■ ■ Xr,\Vf,) , 

(compare (12. 4[) ). The difference, of course, is that this is a canonical ensemble and 
the partition function must be obtained by summing the weights over only those 
configurations with the proper numbers of all species. This relation between this 
truck model on a ring of L + 1 sites and the two species open system studied in this 
paper is completely parallel to that between the usual defect particle model and 
the open one species TASEP. We expect that details of the stationary state could 
be worked out in parallel to that of the usual defect particle model, but we have 
not done so. 

5.2 In the second model, the ring has N sites labeled by i S [—N/2 + l,N/2] 
and contains A^i — piN first class particles, N^ — P2N second class particles, and 
Nq = N — Ni — N2 = PqN holes. The particles jump clockwise according to the 
TASEP rules given in section [H except at one specified semi-permeable bond, say 
between site and site 1, which prohibits the passage of second class particles. (We 
can think of this bond as a restriction in a channel) . 

Unfortunately we do not have an exact solution for this system. To see what 
happens, however, we note that, as in the open system, the current J2 of second 
class particles will vanish in the stationary state. On the other hand, since we 
would have J2 = p2(po — Pi) if the system were uniform, a uniform state is possible 
only if pi — po- If Pi < Po then J2 would be positive in the uniform system and 
second class particles would drift to the right; the upshot is that there will be a fat 
shock of width w = p2N/{l — 2pi) containing all second class particles at density 
P2 = 1 — 2/0i pinned to the back of the barrier. If pi > p2 then the fat shock of 
width w — p2N/{\ — 2po) will be pinned to the front of the barrier. In the case 
Po = Pi = (1 — P2)/2 the system will be uniform. See Figure [H for some typical 
profiles in this system; note that A'o, A^i, and N2 have been chosen so that the bulk 
densities in Figures (a) and (b) here are the same as those in Figures (a) and (e) 
of Figure [21 but that the boundary effects and finite density shock transition are 
noticeably different in the two models. 
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ONE SEMI-PERMEABLE BOND: N„ = 450; N, = 300; N, = 250 ONE SEMI-PERMEABLE BOND: N„ = 375; N, = 375; N, = 250 
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Figure 4. Density profiles in a system with one semi-permeable 
bond: L = 1000. 



Letting iV — > oo with p\ , po fixed would yield an infinite system with a barrier 
at the bond (0, 1). Consider first the case p\ < po- To the right of the origin there 
would be no second class particles and a uniform density of first class particles 
described by a product measure. Far to the left of the barrier the state would 
be i/Po.po^ [q^ g^ uniform state with pi — po — pi and P2 — ^ — 2pi. We do not 
know, however, the structure of the system just to the left of the barrier. Similar 
conclusions hold for pi > po- 
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Appendix A. A particular representation 



A representation of the algebra (|2.2P " (|2.3p which satisfies (|2.6[) may be obtained 
from [H and [9]: 



(A.l) Ai 



/ 1 1 

110 

11 

1 

V ■ 



\ 



Xo 



/ 1 
110 
Olio 
11 

V ■ 



\ 



(A.2) 



X2 — AiAo — AqAi — [Ai, Ao] — 



/ 1 







V ■ 
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(A.3) (VK^I = 1, 



\-a 



a 



1-a 



Wb 



( 1 \ 



1-/3 



V 



/ 



The exponential growth of the components of (Wc\ and |Va) for certain values of 
a and /3 in fact causes no concern here: because we always have n > 0, the matrix 
product needed to calculate the probability of any configuration r (see (12 .4^ ') will 
contain at least one factor X2, and using (|2.6p one can see that this implies that 
the corresponding matrix element is finite. 



Appendix B. Asymptotics of the partition function 

We summarize here the asymptotics of the partition function which are needed 
in Section[5l For the case with no second class particles ^7j we need Z"''^ only when 
a = 1 and/or /3 = 1: 



(B.l) 



Z"'i(j,0)-Zi^"(j,0) 



l-2a / 1 
(l-a)2 \a[\-a) 
2 4^ 



/TTJ 



1/2' 



A? 



+1 



I V^(2a - 1)2 j3/2 ' 



if a< 1/2, 
if a = 1/2, 
if a > 1/2. 



The generating function is 

00 
(B.2) Y.^^'l^- 



.0 



2a 



2/3 



^^^ ^2a-l + x/l-4A/ V2/3-1 + V1-4A 

For the model with second class particles [2]: 
• In region I, (pic < a, /3) 

2 \ i+i 



(B.3) 



na(3\/ L^ — "n? 



AL' 



(B.4) 



. 7rL((2a - l)i + n)((2/3 - l)i + n) \L^ -n' 
• In region II, (a < ac, a < P) 

/3(l-2a) / 1 \^+V a 



L — n 
L + n 



Z''''^{L,n)^ 



a{P — a) \a{l — a) J \1 — a 
• On the boundary of regions I and II, (ac ~ a < (3) 



(B.5) Z"'^(L,n) 



j3n{L — n) 



4L2 



L+l 



L — n 
2L{{2l3-l)L + n) \L^ -n^) [l + n 

• On the boundary of regions II and III, {a — P < ac) 



(B.6) Z"''^(L,n) 



{l-2a){{l~2a)L-n) 
(l-a)2 



1 



Qfl 



a(l ~ a) J (1 — a) 



L+1 



^p^^ ^^_ n{L-n) JL^-n^ f AL^ y^(L-n 
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• At the triple point, {uc = a = j3) 

(^■^) ^"'(^'")- 2L(L + n)V L. KL^-nV \L + n 

Asymptotics in Region III and on the I/III boundary are obtained from those of 
Region II and the I/II boundary by exchange of a and /3. 



Appendix C. Finite volume corrections to density profiles 

We consider here again the problem of finding asymptotic values of the density 
profiles, beginning with a discussion of the method of [2\. The partition function 
can be expressed as 

(C.l) Z'^'''(L, n) = ^ [i?(L, n, P) - R{L, n, a)] , 

a ~ p 

where 

with C™ the Catalan triangle numbers (|3.8p . An asymptotic analysis of (jC.2p then 
leads, through (jC.ip and the formulas p.l3p - (|3.14|) for the densities, to the density 
asymptotics. In 2J the asymptotic density at position x was calculated as 

(C.3) lim {T]aiiL)) ^^.f , , 

with iL — [xL\ . As observed in Section [5l however, if x is the location of the fixed 
shock in regions II or III, and one considers limits as in (|C.3p with i^ = [xL\ +c\/L, 
then the limiting density value depends on c. This c dependence may be calculated 
by the methods of Section [5] (see for example (|5.22p ): here we sketch briefly an 
alternate and more direct method which extends the work of [2]. 

The key step is the computation of the asymptotics of R{L, n, a); it is convenient 
to introduce ac = {L — n)/{2L) (see (|1.2p ). We must determine which terms in 
(|C.2p dominate the sum. If we let L — > oo at fixed n and a there are three regimes: 
(i) a > ac, for which the maximum of the summand is attained when k is of order 
L and the sum can be approximated by a Gaussian integral; (ii) a < ac, in which 
the maximum is attained when k is of order —L and the sum can be approximated 
by a geometric series; and (iii) a = ac, for which the maximum occurs when k is of 
order 1 and the sum can be approximated by half of a Gaussian integral. However, 
there are intermediate regimes in which the sum is dominated by terms in which k 
is of order ± vL, and it is these which generate the finite volume density corrections 
that we seek. 
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One needs an asymptotic estimate of R{L, n, a) which holds for all large L and 
n. Such an estimate is R{L, n, a) ^ RiL, n, a), where 



(C.4) R{L,n,a) 



a(l-a)/ \l-a) \ ^ a(L + n) 

a < ac, 
2nL / 4:L^ 



(n - L{1 - 2a))^7rL(L2 - n^) 



L-nY L(l-2a)-n\ 



y/a{L + n) 

here $ is as 1 in (|5.18p and '^{t) = \/2TTe*^/^\t\^{t). The asymptotic estimate 
R ^ R holds in the sense that for a and the ratio n/L uniformly bounded away 
from and 1 the quantity \R/R — 1| is small when L is large — more precisely, for 
any e > there is a constant C^ such that \R/R— 1| < CeL^^/^^'. We remark that 
the two forms in (IC.4|) in fact agree for ac < a < ac + 0{l/y/L). 

From (jC.l|l and (jC.4p one obtains similarly improved asymptotics for the par- 
tition function Z"'^{L,n), and the full density asymptotics then follows from the 
exact formulas of [^ or Theorem [ 
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